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Abstract 

We extend our calculation of the atmospheric neutrino fluxes to polar and tropical regions. It 
is well known that the air density profile in the polar and the tropical regions are different from 
the mid-latitude region. Also there are large seasonal variations in the polar region. In this 
extension, we use the NRLMSISE-00 global atmospheric model [l| replacing the US-standard ’76 
atmospheric model Q], which has no positional or seasonal variations. With the NRLMSISE-00 
atmospheric model, we study the atmospheric neutrino flux at the polar and tropical regions with 
seasonal variations. The geomagnetic model IGRF Q] we have used in our calculations seems 
accurate enough in the polar regions also. However, the polar and the equatorial regions are the 
two extremes in the IGRF model, and the magnetic field configurations are largely different to 
each other. Note, the equatorial region is also the tropical region generally. We study the effect of 
the geomagnetic field on the atmospheric neutrino flux in these extreme regions. 
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I. INTRODUCTION 


Ml 


In this paper, we extend the calculation of the atmospheric neutrino flux 14- 6| to the sites 

n Q 

in polar and tropical regions. In our earliest full 3D-calculation J^, we used DPMJET-III [7j 
for the hadronic interaction model above 5 GeV, and NUCRIN 8| below 5 GeV. We modihed 


DPMJET-III as in Ref. 


muon spectra better, mainly using 
q, we introduced JAM interaction 


to reproduce the experimenta 
the data observed by BESS group j^. In a recent work 
model for the low energy hadronic interactions. JAM is a nuclear interaction model devel- 

n n 

oped with PHITS (Particle and Heavy-Ion Transport code System) [ij . In Ref. [6| , we could 
reproduce the observed muon flux at the low energies at balloon altitude [m| with DPMJET- 
III above 32 GeV and JAM below better than with the combination of DPMJET-III above 
5 GeV and NUGRIN below that. Besides the interaction model, we have also improved the 
calculation scheme according to the increase of available computational power, such as the 


■‘virtual detector correction” introduced in Ref. [5[ and the optimization of it in Ref. jb]. 
The statistics of the Monte Garlo simulation is also improved at every step of the work. 

We used the US-standard ’76 atmospheric model in our earlier calculations. The US- 
standard ’76 atmospheric model had been used generally in the study of cosmic rays in the 
atmosphere for a long time 13|. But, the air density profile in US-standard ’76 is repre¬ 
sented as a function of altitude only, and has no time variation and no position dependence 
around the Earth. In Ref |l^, we discussed the validity of using such a atmospheric model in 
the calculation of atmospheric neutrino flux, assuming small variations to the US-standard 
’76 atmospheric model. However, the difference of air density prohle in the polar regions 
and mid-latitude regions is lager than the considered variations in the study. Also there 
is a large seasonal variation of the air density prohle in the polar region. We install the 


NRLMSISE-00 global atmospheric model [l| which represents proper position dependence 
and the time variations on the Earth, to calculate the atmospheric neutrino hux in the polar 
and tropical regions. 

We have used the IGRF geomagnetic held model in our calculation, and it is accurate 
enough in the polar and tropical (~ equatorial) regions. However, the geomagnetic held 
strongly ahects the atmospheric neutrino hux, and is largely diherent in the polar and 
equatorial regions. The extension in this paper is also the study of atmospheric neutrino 
hux under these widely diherent geomagnetic held conditions. 
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The models of primary cosmic rays spectra and interactions are the same as those in 


Ref. 

BESS 


GL The model of primary cosmic rays are constructed based on the AMSOl [1^ and 
161 observations. However, there are newer cosmic ray observation experiments [ 3 - 


3| . and we will make a short comment on them and the error of our calculation in the 
summary. We note, the combination of the modihed DPMJET-III above 32 GeV and JAM 
below that reproduces the observed muon spectra best with the present cosmic ray spectra 
model. 


In our 3D-calculations of the atmospheric neutrino flux, we followed the motion of all the 
cosmic rays, which penetrate the rigidity cutoff, and their secondaries. Then we examine 
all the neutrinos produced during their propagation in the atmosphere, and register the 
neutrinos which hit the virtual detector assumed around the target neutrino observation site. 
Therefore, we do not need a change in the calculation scheme other than the atmospheric 
model to calculate the atmospheric neutrino flux at a new site in the polar and tropical 
regions. 

We study in detail the atmospheric neutrino flux at the India-based Neutrino Observatory 
(INO) site (lat, lon) = (9°59", 77°16") for the tropical and equatorial region, and the South 
Pole (—90°00", 0°00") and Pyhasalmi (63°40", 6°41") mine (Finland) for the polar regions in 
this paper. Also we compare the atmospheric neutrino flux calculated with the NRLMSISE- 
00 atmospheric model and that calculated with the US-standard ’76 atmospheric model at 
Super Kamiokande (SK) site (36°26", 137°10"). 

We note that the atmospheric neutrino production height is important for the anal¬ 
ysis of neutrino oscillations as well as the flux. As the production height is related to 
the atmospheric model, we study the production height of atmospheric neutrinos with the 
NRLMSISE-00 atmospheric model in detail, and compare the atmospheric neutrino pro¬ 
duction height with that calculated with US-standard ’76 atmospheric model at the SK 
site. 


The tables for the flux and production height calculated in this paper are available in the 
web page, http://www.icrr.u-tokyo.ac.jp/~mhonda . 
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II. NRLMSISE-00 ATMOSPHERIC MODEL 


NRLMSISE-00 [l| is an empirical, global model of the Earth’s atmosphere from ground to 
space. It models the temperatures and densities of the atmosphere’s components. However, 
the air density prohle is the most important quantity in the calculation of atmospheric 
neutrino flux. We calculate the ratio of the air density in 4 seasons, March - May, June - 
August, September - November, and December - February at the SK site (KAM), INO site 
(INO), South pole (SPL), and Pyhasalmi mine (PYH) by the NRLMSISE-00 atmospheric 
model to that by the US-standard ’76 atmospheric model, and show it in Fig. [T] as a function 
of altitude. 




Altitude (km) 



Altitude (km) 



Attitude (km) 


Attitude (km) 


FIG. 1. The ratio of air density in the NRLMSISE-00 atmospheric model to that of US-standard 
’76 atmospheric model for the SK site (KAM), INO site (INO), South pole (SPL), and Pyhasalmi 
mine (PYH), in the 4 seasons, March - May, June - August, September - November, and December 
- February. 


In the tropical region (INO), we see a ~ 20 % larger air density than the US-standard ’76 
atmospheric model at the altitude of ~ 15 km a.s.l (above sea level). However, except for 
that, the air density prohle is similar to the US-standard ’76 atmospheric model, and there is 
almost no seasonal variation by the NRLMSISE-00 atmospheric model. In the mid-latitude 
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region (KAM), we find seasonal variations but they are small, and the air density prohle is 
close to that of the US-standard ’76 atmospheric model through all the seasons below 40 km 
a.s.l. On the other hand, we find large seasonal variations in the Polar region (SPL and 
PYH) above 10 km a.s.l., especially at the South Pole (SPL). At the Pyhasalmi mine, the 
air density prohle is similar to that of the US-standard ’76 atmospheric model below 10 km 
a.s.l, but at the South Pole, the air density decreases quicker than that even below 10 km 
a.s.l. 

Thus, we expect some seasonal variations of atmospheric neutrino hux except for the 
tropical region, and study it some in detail with the NRLMSISE-00 atmospheric model in 
the section El 


III. GEOMAGNETIC FIELD AND SITES 


We have already reported a large effect of geomagnetic held on the calculation of a ;- 
mospheric neutrino hux for several sites in mid-latitude, such as SK and SNO sites Q, 6|, 
through the rigidity cutoh and muon bending. 

As a naive illustration of the rigidity cutoh, consider the gyro motion of cosmic rays 
guided by the horizontal component of the geomagnetic held near the Earth. The cosmic 
rays with a very small radius of gyro motion can not arrive at a very close point to the earth. 
As the radius of the gyro motion becomes larger, the cosmic rays can access the point, but 
the Earth works as a slant shield, limiting the access azimuth angle of cosmic rays. When 
the radius is large enough, the Earth becomes a hat shield just limiting the upward going 
cosmic rays, and the limitation on the access azimuth angle disappears. This mechanism is 
called the rigidity cutoh. 

The ehect of the muon bending on the atmospheric neutrino hux is often explained by the 
diherence of the arrival zenith angle of neutrinos from that without the muon bending by 
the geomagnetic held. As the atmospheric neutrino hux above a few GeV has a large arrival 
zenith angle dependence, a little diherence of the zenith angle results in a large diherence 
of the hux separately visible from the rigidity cutoh at these energies. For the diherence of 
the arrival zenith angle, the horizontal component of the geomagnetic held is responsible. 

Thus, the horizontal component of geomagnetic held (i?h) is an important parameter to 
understand the ehects of rigidity cutoh and muon bending. In Fig. [2l we draw the strength 
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of the horizontal component of the geomagnetic held using the IGRF geomagnetic model for 
the year of 2010 with the position of the SK site {Bh ~ SOOOOnT), INO site {Bh ~ 40000nT), 
South Pole {Bh ~ lOOOOnT), and Pyhasalmi mine {Bh ~ ISOOOnT) for which we are going 
to calculate the atmospheric neutrino hux. 


■ SK ■ INO 4 Pyhasalmi 



FIG. 2. The horizontal component and the sites, Kamioka (KAM), South pole (SPL), INO site 
(INO), and Pyhasalmi mine (PYH). where atmospheric neutrino flux is calculated. 


IV. CALCULATION OF ATMOSPHERIC NEUTRINO FLUX 


Except for the atmospheric model, the calculation scheme, including the interaction model 
and the cosmic ray spectra model are the same as in the previous work 6|. We assume the 
surface of the Earth as a sphere with radius of = 6378.140 km. In addition, we assume 
2 more spheres, the injection sphere with a radius of Rinj = Re + 100 km and escape sphere 
with radius Resc = 10 x Rg. Before Ref. j^, we assumed one more, the simulation sphere 
with radius of Rsim {Rinj < Rsim < Resc)- We discarded the cosmic rays which go outside of 


this sphere after the injection. We took Rsim = Re R 3000 km in Ref. However, now we 
identify the simulation sphere and the escape sphere by taking Rsim = Resc after Ref. jo]. 

For each cosmic ray event simulation, we sample an energy and a chemical composition 
of cosmic ray to simulate, according to the cosmic ray spectra model. Then we sample 
the position and the initial direction of the cosmic ray on the injection sphere to start the 
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simulation. For each cosmic ray, we apply the rigidity cutoff test and check if it reaches this 
position penetrating the rigidity cutoff. 

This rigidity cutoff test is carried out by the back tracing method, solving the equation 
of motion in the inverse time direction in the geomagnetic held. When the sampled cosmic 
ray reaches the escape sphere without touching the injection sphere again, we judge that the 
cosmic ray can reach the starting position from deep space, and feed it into the simulator of 
the propagation in the atmosphere. 

Generally, the transition from inhibited to allowed rigidity is not clear. Most cosmic rays, 
which fail the rigidity cutoff test, hit the injection sphere very quickly before completing one 
cycle of gyro motion. However, some of them with rigidity near the transition, travel a long 
distance before hitting the injection sphere. Also a cosmic ray with slightly higher rigidity 
than one which passes the rigidity cutoff test may fail the test even starting from the same 
position and direction. 

A primary cosmic ray may pass through the injection sphere more than once. In this case 
we have to select the first entrance to avoid double counting. For this purpose, we require 
that the cosmic ray should not come back to the injection sphere again during back tracing. 

The simulator of cosmic ray propagation in the atmosphere follows the motion of all the 
cosmic rays, both primary cosmic rays which passed the rigidity cutoff test and secondary 
cosmic rays produced in interactions of the cosmic rays recursively, until they reach the 
escape sphere, or hit the surface of the Earth, or interact with an air nucleus, or decay. 
Each neutrino produced in the cosmic ray interaction or decay is examined if its path will 
take it through the virtual detector assumed around the neutrino observation site, and when 
it goes through the virtual detector, it is registered and the number is used to calculate the 
atmospheric neutrino flux. We take a circle with radius of 1113.2 km as the virtual detector. 
Note, the radius corresponds a change in longitude of 10 degrees on the equator. 

As the virtual detector is far larger than the real neutrino detector, we introduce the 
“virtual detector correction”. Assuming a circle with a radius of 9^ around the real detector 
as the virtual detector, the flux dehned as the average flux in the virtual detector may 
be written in the form as 

•hd = $0 + ^0 + • • •, (1) 

with the flux $0 at the real detector j^. Then, we can cancel out the term, using two 


fluxes $1 and $2 determined in the virtual detectors with radii 61 and 62 respectively as 



~ el-el ~ 1 -r2 


( 2 ) 

[g] to 


, but optimize it to r = l/\/2 in Ref. 


minimize the statistical error. 

V. ATMOSPHERIC NEUTRINO FLUX AT EACH SITE 

In Fig. [3l we show the one year average of atmospheric neutrino fluxes at the SK site, 
INO site, South Pole, and Pyhasalmi mine, averaging over all the directions. The qualitative 
features are the same at all the sites, but we find a difference of flux among the sites by 
factor ~ 3 at the low energy end due to the large difference of the cutoff rigidity among these 
sites. The differences of flux among the sites above 10 GeV is small in the figure. However 
we expect a large seasonal variations at South Pole and Pyhasalmi mine, and we will study 
them in the next subsection some in detail. 

A. Seasonal variation of the atmospheric neutrino flux 

To study the seasonal variations, we calculate the ratio of seasonal fluxes to the yearly 
flux average, and show them in Fig. 01 Also the ratio of the flux calculated with the US- 
standard ’76 atmospheric model to the yearly flux average is shown in the panel for SK site 
below 1 TeV. Note, the calculation of the fluxes with the US-standard ’76 atmospheric model 
was carried out in the previous work , and the statistics of the Monte Carlo simulation is 
poorer than that in the present work especially above 1 TeV. 

At the mid-latitude region (SK site) and tropical region (INO site), the seasonal differ¬ 
ences are small, and is difficult to see even in the ratio. On the other hand, the seasonal 
variaston is large in the Polar region, as the summer-winter difference reaches more than 
10 % at the South Pole at 10 TeV, and more than 5 % at the Pyhasalmi mine. At both sites, 
the atmospheric neutrino flux is higher in summer (December - February at the South Pole 
and June - August at the Pyhasalmi mine). This may be understood by the fact that the 
air density at higher altitudes (> 15 km) is higher in the summer at both sites (see Fig. [T]). 
We note the seasonal variation of the 1 /^ and i>e fluxes starts from < 10 GeV, where that of 
the and is still small. 
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The effect of the air density prohle on the atmospheric neutrinos is generally discussed 
considering the relative probability of pions to decay or to interact with air nuclei. This effect 
works equally on all flavors of neutrino at high energies (> 100 GeV), where the probability 
of the parent pions to interact becomes comparable to that to decay, and may explain the 
large seasonal variation above 1 TeV. However, it does not explain the seasonal variation of 
the Ve and v^. fluxes starting from < 10 GeV, and the difference from that of the and 
fluxes. 

To explain these difference, we need to consider the muon propagation in the atmosphere. 
When the air shrinks down lower, the muons are created at lower altitudes, and the proba¬ 
bility of muons to hit the ground before decaying increases. When the muons hit the rock or 
ice, they lose their energy quickly producing neutrinos with very low energies (< 0.1 GeV) 
only. Then the flux of neutrinos produced in the decay of muons decreases. This mechanism 
results in the variation of neutrino flux near the vertical directions at relatively high energies 
(> 10 GeV). It explain the variation of neutrino flux there, and the difference of the Vf. and 
Vf. fluxes and the and fluxes, as half of the and z/^ are created in the pion decay 
directly. 

There is another mechanism which has a large effect on the energy spectra of neutrinos 
at lower energies. When the atmosphere shrinks lower, the muons travel in denser air and 
lose more energy before the decay. The larger energy loss of the muons causes a shift to 
the energy spectra of the neutrinos produced in the muon decay towards the lower energy 
direction. As the energy spectra of atmospheric neutrinos are steeply decreasing at the 
energies ^0.1 GeV, the fluxes decrease in the denser air at those energies. However, this 
mechanism is effective when most muons decay in the air (E^ < GeV), and for the neutrinos 
with energies below a few GeV. This mechanism is responsible to the seasonal variation of 
neutrino flux below 10 GeV. 


The differences between the fluxes calculated with the US-standard ’76 atmospheric model 
and yearly averaged one at the SK site are far less than 5 %. It is difficult to say that 
the differences are due to the difference in the atmospheric model, since there are some 
improvement in the calculation scheme and statistics of Monte Garlo simulation after the 


previous work 


M- 
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FIG. 3. All-direction averaged atmospheric neutrino flux for 4 sites averaging over one year. KAM 
stands for the SK site, INO for the INO site, SPL for the South Pole, and PYH for the Pyhasalmi 
mine. 

B. Flavor-ratio of the atmospheric neutrino flux 

In Fig. ISl we show the flavor-ratio dehned as the flux ratio of different neutrino flavors 
such as, (z/^ -|- i)^)/{ve + he), and (r'e/he) at the SK site, INO site. South Pole, 

and Pyhasalmi mine, averaging over all the directions. We And that the flavor-ratio is very 
similar to each other among these sites, conhrming the stability of the flavor-ratio. However, 
the flavor-ratio is an important quantity in the study of neutrino oscillations, so we need to 
study the seasonal variations and position dependence’s more precisely. 

To see the seasonal variation of the flavor-ratio, we calculate the ratio of the flavor-ratio 
calculated with the seasonally averaged fluxes to that calculated with the yearly averaged 
fluxes and plot them in Fig. [6] for each site. Also, to see the positional dependence on the 
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FIG. 4. Ratios of the all-direction averaged flux in June - Augnst and in December - February 
to that of the yearly average. For the SK site, we also plot the ratio for the calculation with the 
US-standard ’76 atmospheric model to the yearly average in dash-dot below 1 TeV. 

earth, we take the flavor-ratio calculated with the yearly averaged fluxes at the SK site as 
the “reference flavor-ratio”. Then we calculate the ratio of the “reference flavor-ratio” to 
those calculated with yearly averaged fluxes at other sites, and plotted them in their panels 
of Fig. [6l While in the panel of the SK site, we show the ratio of the flavor-ratio calculated 
with the US-standard ’76 atmospheric model to that calculated with the yearly averaged 
fluxes at the SK site (he. the “reference flavor-ratio”). 

At the sites in the Polar region (South Pole and Pyhasalmi mine), the flavor-ratio, 
(z/^ -h h^)/(r'e + he) shows a seasonal variation, high in summer and low in winter, with 
the maximum of the amplitude at ~ 100 GeV. This is considered to be due to the seasonal 
variation of the altitude of cosmic ray interactions. Also the flavor-ratio shows some differ¬ 
ences from the SK site at the energy below a few GeV. This is considered to be due to the 
lower air density at the neutrino production height of 10~20 km a.s.l (see section IVTl) in the 
Polar region. The smaller muon energy loss causes a smaller shifts of the energy spectra of 
the neutrinos produced in the muon decay . 

Note, the air density at the South pole at the 10~20 km a.s.l. is much lower than at the 
Pyhasalmi mine in Fig. [H but the difference of the flavor-ratio, (z/^ -|- z/^)/(z/e + i^e) at the 
South Pole from the SK site is similar to that at the Pyhasalmi mine. This is considered to 
be due to the higher observation site at the South Pole (2835m a.s.l.). The shorter distance 
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from the ground to the production height of muons reduces the neutrino fluxes produced by 
muon decay at the South Pole. 

In the Ve/^e ratio, we also hud a difference from the SK site at the South Pole and 
Pyhasalmi mine, below a few GeV. This difference is considered to be due to the difference 
of cutoff rigidity. The Ve/^e ratio reflects the /'k~ ratio of parent pions. As the majority 
of primary cosmic rays are protons, there is a 7r+ excess generally. Especially when the 
cutoff rigidity is low enough, the pion production of primary cosmic rays overwhelms that of 
secondary cosmic rays, and and ratios are high even at low energies. However, 

when cutoff rigidity is high, the pion production by secondary cosmic rays can not be ignored, 
and the tt’*' excess is diluted by the secondary neutron cosmic ray interactions. 

In the comparison of neutrino flavor-ratio between the SK site (mid-latitude region) and 
the INO site (tropical region), we hnd a small difference in the + h^)/(r^e + ^e) ratio 
due to the difference of air density at 15 km a.s.l., and the difference of the muon energy 
loss. Other ratios are quite similar to each other. The differences of the flavor-ratio with 
the US-standard ’76 atmosphere model are very small to that in the present calculation. 



FIG. 5. Neutrino flavor-ratio calculated with the all-direction and one-year averaged atmospheric 
neutrino flux. KAM stands for the SK site, INO for the INO site, SPL for the South Pole, and 
PYH for the Pyhasalmi mine. 
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FIG. 6. Ratios of the flavor-ratio calculated with the all-direction and seasonally averaged flux 
to that calculated with all-direction and yealy averaged flux. The solid lines are the flux ratio in 
June - August, and dashed lines are in December - February. KAM stands for the SK site, INO 
for the INO site, SPL for the South Pole, and PYH for the Pyhasalmi mine. Taking the yealy 
averaged flavor-ratio at SK-site as the ’‘reference flavor-ratio”, the ratio of ’‘reference flavor-ratio” 
to the yearly averaged flavor-ratio at each site are plotted with dash-dot other than the panel for 
the SK-site. In the panel for the SK-site, the ratio of the flavor-ratio with the US-standard ’76 
atmospheric model to the yearly averaged flavor-ratio at SK-site is plotted with dash-dot. 


C. Zenith angle variation of the atmospheric neutrino flux at 3.2 GeV 


In Fig. 13 we show the the zenith angle dependence of the atmospheric neutrino flux 
calculated at the 4 sites at 3.2 GeV. We also plot the atmospheric neutrino flux calculated 
with the US-standard ’76 atmospheric model 6(] at the SK site, although the differences are 
small and within the thickness of the line. The seasonal variation is very small at the SK 
and INO sites as in the all-direction average, but is seen clearly in the down going neutrino 
fluxes at the South Pole and Pyhasalmi mine even at this energy. At both the sites the trend 
is the same as the all-direction average at high energies, the atmospheric neutrino flux is 
large in summer (June - August at Pyhasalmi mine and December ~ February at the South 
Pole). 

The amplitude of the seasonal variation is different among different neutrino flavors in 
Fig. 0 It is ~ 2.5% for and z/^, and ~ 10 % for z/g and z/g for the vertically down going 
direction at the South Pole. The muon energy loss in the atmosphere before the decay seems 
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to be mainly responsible to this seasonal variations. 


' " ' 1 " ' 

_PYH 

. 3.2 GeV A 

' 1 "'' 1 " " 




\ V® 

/ / 

ve\' 


Jun-Aug 

Dec-Feb 

, 1 ,,,, 1 ,,,, 




FIG. 7. Arrival zenith angle dependence of atmospheric neutrino fluxes averaged over all azimuth 
angles at 3.2 GeV. KAM stands for the SK site, INO for the INO site, SPL for the South Pole, and 
PYH for the Pyhasalmi mine. The solid lines are the time average in June - August, and dashed 
lines are in December - February. For SK site, we also plot the calculation with the US-standard 
’76 atmospheric model using a dotted line (but is difficult to see separately). Note, 6^ stands for 
the zenith angle. 

Leaving from the seasonal variation, the atmospheric neutrino flux at near horizontal 
directions are largely different from site to site. The horizontal fluxes at the SK site are 
~ 10 % smaller than those at the South Pole or Pyhasalmi mine, but ~ 10 % larger than 
those at the INO site. The differences in the atmospheric neutrino fluxes at near horizontal 
directions is mainly caused by the difference of rigidity cutoff at near horizontal directions. 
The rigidity cutoff works most strongly for the neutrinos arriving from near horizontal 
geomagnetic East directions. At the SK site, we can almost ignore the rigidity cutoff at near 
vertical directions, but it is still strong at near horizontal directions. We can explain the 
difference between the SK site and the Polar region by the difference of rigidity cutoff, but 
can not explain the difference between INO and SK sites only by the difference of rigidity 
cutoff. We need to consider the muon bending, as the suppression for the and z/g is 
stronger than that for and Ve at the INO site. The suppression by the muon bending 
will be understood better with the azimuthal variation of atmospheric neutrinos described 
in subsection IV D[ 
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The fluxes at the South Pole and Pyhasalmi mine are a little smaller than the SK and 
INO sites at cos 6^2 ~ —0.7, and they have some structure there due to the rigidity cutoff 
and the muon bending at the places where the neutrinos are produced. On the other hand, 
the rigidity cutoff and the muon bending effects are weaker at the SK and INO sites at 
cosOz ~ —0.7 than the South Pole and Pyhasalmi mine. 

Apart from the geomagnetic held, the vertically down going neutrino fluxes at the South 
Pole are a little smaller than those at the Pyhasalmi mine due to the difference of the 
observation altitudes. For the South Pole, 2835m a.s.l. is assumed as the observation site, 
and sea level for Pyhasalmi mine. 

D. Azimuthal variation of atmospheric neutrino flux 






FIG. 8. Arrival azimuth angle dependence of the time averaged atmospheric neutrino flux in 5 
zenith angle bins, 1 > cos 02 > 0.6, 0.4 > cos 02 > 0.2, 0.2 > cos 02 > —0.2, —0.2 > cos 02 > —0.6, 
and —0.6 > cos 02 > — 1 at the SK site. The azimuth angle is measured counter-clockwise from 
South. Solid line shows the averaged fluxe in 1 > cos 02 > 0.6, dashed line that in 0.4 > cos 02 > 0.2, 
dotted line that in 0.2 > cos 9z > —0.2, dash dot that in —0.2 > cos 9z > —0.6, and dash 2dots aht 
in —0.6 > cos 02 > —1. 

Next, we study the azimuthal variation of the neutrino fluxes in 5 zenith angle bins; 
1 > cos 02 > 0.6, 0.4 > cos 02 > 0.2, 0.2 > cos 02 > —0.2, —0.2 > cos 02 > —0.6, and 
—0.6 > cos 02 > —1 at 3.2 GeV, averaging over a year. In Fig. [8|, we show the azimuthal 
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variation of atmospheric neutrino flux at the SK site, in Fig. |9l those at the INO site, in 
Fig. [TOl those at the South Pole and in Fig. [TTl those at the Pyhasalmi mine. 

In these hgures we observe two kinds of effects from the geomagnetic held. One is the 
rigidity cutoff, and the other is the muon bending. The geomagnetic held directed to the 
North hlters the low energy cosmic rays from the East directions, since the cosmic rays 
generally carry a positive charge. The rigidity cutoh reduces all havors of neutrino from the 
East at the same rate. 

On the other hand, the ehect of muon bending depends on the muon charge. The positive 
muon and its decay products are ahected by the geomagnetic held in the same way as the 
rigidity cutoh, and the negative muon and its decay products are ahected in the opposite 
way. Then, muon bending reduces the and huxes from the East, but reduces the 
and he huxes from the West. 

In Fig. [HI we hnd small azimuthal variations of atmospheric neutrino huxes at near 
vertically down going directions (1 > cos 6^2 > 0.6) at the SK site. As the variation shape 
is almost the same among all neutrino havors, this variation is considered to be due to the 
rigidity cutoh, but the ehect is already small at the SK site for near vertical directions at 
3.2 GeV. 

At near horizontal directions ( 0.2 > cos^^ > ~ 0 . 2 ), and Vf. huxes have large sinusoidal 
azimuthal variations, but and z/g huxes have small but more complicated azimuthal varia¬ 
tions. This is because the rigidity cutoh and the muon bending work in the same directions 
for and v^. huxes, and in the opposite directions for and z/g huxes. The dip at ~ 90° 
commonly seen in all neutrino havors is considered to be due to the rigidity cutoh, and the 
dip at ~ 270° seen in the and z^e huxes is due to muon bending. 

In Fig. ini we hnd larger azimuthal variations of atmospheric neutrino huxes at near 
vertically down going directions (1 > cos 6*2 > 0.6) at the INO site compared to the SK site. 
As the variation shape is similar among all neutrino havors, this is considered as an ehect 
of the rigidity cutoh, and is still strong at the INO site for vertical down going directions 
at 3.2 GeV. The azimuthal variation at near horizontal directions (0.2 > cos 6^2 > —0.2) are 
also larger than the SK site for each neutrino havor, due to the larger horizontal component 
of the geomagnetic held. 

We note that the muon bending ehect is seen in the azimuthal averaged hux plot (Fig. [7]) 
at the INO site at near horizontal directions. The muon bending suppresses the and v^. 
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FIG. 9. Azimuth angle dependence of the time averaged atmospheric neutrino flux in 5 zenith 
angle bins, 1 > cos 02 > 0.6, 0.4 > cos 02 > 0.2, 0.2 > cos 02 > —0.2, —0.2 > cos 02 > —0.6, and 
—0.6 > cos 02 > —1 at the INO site. The azimuth angle is measured counter-clockwise from the 
South. Solid line shows the averaged fluxe in 1 > cos 02 > 0.6, dashed line that in 0.4 > cos 02 > 0.2, 
dotted line that in 0.2 > cos 02 > —0.2, dash dot that in —0.2 > cos 02 > —0.6, and dash 2dots aht 
in —0.6 > cos 02 > —1. 

fluxes from the West, but enhances them from the East. However, the rigidity cutoff works 
strongly for the East directions, canceling the enhancement from the East, and the muon 
bending is seen as a suppression of the and he fluxes in Fig. [3 On the other hand for 
and Ue, the muon bending enhances the fluxes from the West, but suppresses them from the 
East. As the rigidity cutoff works weakly for the West directions, the muon bending is seen 
as an enhancement of the and z/g fluxes in Fig. [3 The same mechanism works at other 
sites, but the amplitudes of the suppression and enhancement are small, and it is not seen 
clearly even at the SK site (KAM). 

In Fig. [ini we hnd that there are almost no azimuthal variations of the atmospheric 
neutrino flux except for the near horizontal directions (0.2 > cos 02 > —0.2) at the South 
Pole, since the geomagnetic held is almost vertical at the South Pole. At near horizontal 
directions we hnd the azimuthal variations of the and having the minimum at ~ 90° 
and the maximum at ~ 270°. Also there are slightly smaller variations of z/^ and z/g, having 
the minimum at ~ 90° and the maximum at ~ 270°. These features represent the fact 
that the rigidity cutoh works very weakly at the South Pole, and the muon bending works 
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FIG. 10. Azimuth angle dependence of the time averaged atmospheric nentrino flnx in 5 zenith 
angle bins, 1 > cos 02 > 0.6, 0.4 > cos 02 > 0.2, 0.2 > cos 02 > —0.2, —0.2 > cos 02 > —0.6, and 
—0.6 > cos 02 > — 1 at the South Pole. The azimuth angle is measured counter-clockwise from the 
meridian line of 180 degree in longitude. Solid line shows the averaged fluxe in 1 > cos 02 > 0.6, 
dashed line that in 0.4 > cos 02 > 0.2, dotted line that in 0.2 > cos 02 > —0.2, dash dot that in 
—0.2 > cos 02 > —0.6, and dash 2dots aht in —0.6 > cos 02 > —1. 


mainly with the residual horizontal component of the geomagnetic held. We note that the 
azimuthal variation for upward going directions (cos 02 < —0.2) is also very small at the 
South pole. This is probably because the South Pole is close to the geomagnetic South pole 
in the dipole approximation of the geomagnetic held. 


In Fig. [m we hnd the features of azimuthal variations of down going and near horizontal 
atmospheric neutrino hux (cos 02 > —0.2) at the Pyhasalmi mine are almost the same as 
those at the South Pole. We note that the horizontal component of the geomagnetic held 
at the Pyhasalmi mine {B^ ~ 13000 nT) is even smaller than that of the South Pole ~ 
16000 nT). However, as the Pyhasalmi mine sits at a distant position from the geomagnetic 
North Pole, the rigidity cutoh and muon bending ahects the azimuthal variations of the 
upward going neutrinos. 
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FIG. 11. Azimuth angle dependence of the time averaged atmospheric nentrino flnx in 5 zenith 
angle bins, 1 > cos 02 > 0.6, 0.4 > cos 02 > 0.2, 0.2 > cos 02 > —0.2, —0.2 > cos 02 > —0.6, and 
—0.6 > cos 02 > — 1 at the Pyhasalmi mine. The azimnth angle is measnred counter-clockwise 
from the Sonth. Solid line shows the averaged fluxe in 1 > cos 02 > 0.6, dashed line that in 
0.4 > cos 02 > 0.2, dotted line that in 0.2 > cos 02 > —0.2, dash dot that in —0.2 > cos 02 > —0.6, 
and dash 2dots aht in —0.6 > cos 02 > —1. 

VI. PRODUCTION HEIGHT OF ATMOSPHERIC NEUTRINO 

Here we study the production heights of atmospheric neutrinos calculated with the 
NRLMSISE-00 atmospheric model. The production height is an important parameter for 
analyzing the neutrino oscillations using atmospheric neutrinos. As the production heights 
are distributed in a wide range from sea level to ~ 100 km a.s.h, we examine the cumulative 
distribution, and plot the height where the cumulative distribution of production height 
reaches 10 %, 50 %, and 90 % as a function of neutrino energy for vertical down going neu¬ 
trinos (cos 02 > 0.9) and horizontally going neutrinos (0.1 > cos 02 > 0) in Figs [121 [131 [151 
and [m In those hgures, the heights shown by the lines of 10 %, 50 %, and 90 % indicate 
10 %, 50 %, and 90 % of neutrinos are created below those height, respectively, for the SK 
site, INO site. South Pole, and Pyhasalmi mine respectively. To make the distributions in 
the hgures, we combine particle and anti-particle neutrinos, and average over all azimuth 
angles. In the hgures for the SK site, we also plot the same havor-ratio with the US-standard 
’76 atmospheric model below 32 GeV, calculated in the previous work jo]. 
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The production heights are largely different between and Ue + Pe, and also between 

vertical and horizontal directions. The large energy dependence of production height of the 
Ug and Pe are explained by the fact that they are mainly produced by muon decay in the 
energies 100 GeV for near vertical directions, and < 1 TeV for near horizontal directions. 
The muons fly long distances before decays, and make neutrinos at a lower altitude. As 
the muon flight length increases with the muon energy, the production height of the z/g and 
Pe becomes lower as long as the muon is their major source. However, above 100 GeV for 
vertical directions and 1 TeV for horizontal directions, Kaon decay becomes the main source, 
then the production height becomes high again. 

Muon decay is also a source of z/^ and P^. A similar variation to Ug and Pe is also seen in 
i>fj_ and z/^, but with a smaller amplitude. The fraction of and z/^ created in muon decay 
is not so large, roughly a half at low energies and decreasing with the energy as muons tend 
to hit the ground. The effect of muon flight in the production height is smaller on z/^ and 

than that on Ug and Pg- 

Note, above 100 GeV for vertical directions and 1 TeV for horizontal directions, Kaon 
decay becomes the main source of all flavors of neutrinos. Then the production height 
becomes high, close to the production height of Kaons, or the hadronic interaction zone of 
cosmic rays. We would also like to note that the production height calculated in 3D-scheme 
is smoothly connected to that calculated in the ID-scheme as well as the flux. 

At the SK and INO sites, we observe small difference between the averages of June - 
August and December - February. Especially at the INO site, the difference is almost 
invisible in the figures. The difference between the calculations with the NRLMSISE-00 
and US-standard ’76 atmospheric models at the SK site is also small. On the other hand, 
we hnd a large seasonal variation of production height at the South Pole. The variation 
of the median (cumulative distribution of 50 %) height extend to ~ 20% variation at near 
horizontal directions, and the variation of the amplitude is similar between z/^ -|- and z/g 
-|- Pg. At the Pyhasalmi mine, we also find a seasonal variation of ~ 10 % for the median. 

The energy dependence of the production height is similar among the different sites. 
However, the of 10 % line of the cumulative distribution for z/g -|- Pg at the South Pole is 
a little higher than that at other sites due to the altitude of the sites, as we assumed the 
observation site is 2835m a.s.l. for the South Pole, and sea level for other sites. 

In Fig. [16] we plotted the azimuthal variation of the median of the production height at 
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FIG. 12. Production height of atmospheric neutrinos for and I'e + going to vertically 

down and horizontal directions, summing all azimuth directions at the SK site. The height that 
cumulative distribution reaches 10 %, 50 %, and 90 % are shown as the function of neutrino energy. 
The time average values in June - August are shown as solid lines and in December - February as 
Dashed line. The same values with the US-standard ’76 atmospheric model are also shown in dash 
dot in the 3D region. 

near horizontal directions for the atmospheric neutrinos at = 3.36 GeV (3.16 < Ey < 
3.55 GeV), summing the production height distribution over a year. We hnd there are 
sinusoidal variations with azimuth angle for all flavor neutrinos, but in an opposite direction 
to each other among the neutrinos and anti-neutrinos, at the SK and INO site. The azimuth 
variation of production height at the Pyhasalmi mine is smaller than that at the SK and 
INO sites, but the shape is similar. The azimuthal variation of production height at the 
South Pole is small. These sinusoidal features at the SK and INO sites could be understood 
if we consider that the production height is mainly controlled by the horizontal component 
of the geomagnetic field and the effect of the rigidity cutoff is small. 


VII. SUMMARY AND DISCUSSION 

We have extended our calculation of atmospheric neutrino flux to the tropical (INO site) 
and polar regions (South Pole and Pyhasalmi mine). In this extension, we have updated the 
atmospheric model from the US-standard ’76 model to the NRLMSISE-00 model. When we 
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FIG. 13. Production height of atmospheric neutrinos for and i^e + going to vertical 

down and horizontal directions, summing all azimuth directions at the INO site. The height that 
cumulative distribution reaches 10 %, 50 %, and 90 % are shown as a function of neutrino energy. 
The time average values in June - August are shown as solid lines and in December - February as 
Dashed line. 


compare the two atmospheric models, we hnd the NRLMSISE-00 atmospheric model almost 
agrees with the US-standard ’76 atmospheric model in the tropical and mid-latitnde region. 
However, they disagree with each other largely in the polar region. Also the NRLMSISE-00 
atmospheric model snggests a large seasonal variation in the polar region. The atmospheric 
neutrino flux calculated with the NRLMSISE-00 model is compared with that calculated 
with the US-standard ’76 model at the SK site. 

Adding to the air density prohle, there is a large difference of the geomagnetic held 
conhguration between equatorial and polar regions. It is well known that the inhuence of 
geomagnetic held on the atmospheric neutrino is large. The extension in this paper is also 
the study of the atmospheric neutrino hux in the two extremes in the IGRF geomagnetic 
held model. Note, our calculation so far was limited to the sites in mid-latitude region. 

As expected, the calculated atmospheric neutrino hux at the equatorial (tropical) site 
(INO) suhers a strong ehect from the horizontal component of the geomagnetic held. A 
large reduction of the neutrino hux is still seen at 3.2 GeV for down going directions due to 
the rigidity cutoh. The muon bending causes a large azimuthal variation of neutrino huxes. 
reducing and huxes and enhancing a little and hues at horizontal directions. 
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FIG. 14. Production height of atmospheric neutrinos for and I'e + going to vertically 

down and horizontal directions, summing all azimuth directions at the South Pole. The height that 
cumulative distribution reaches 10 %, 50 %, and 90 % are shown as a function of neutrino energy. 
The time average values in June - August are shown as solid lines and in December - February as 
Dashed line. 

On the other hand, the large seasonal variation of air density profile at the polar regions 
(South Pole and Pyhasalmi mine) causes a large seasonal variation of atmospheric neutrino 
flux. We had expected that the seasonal dependence is large at higher energies since the vr 
decay and interaction ratio changes with the air density. But at the South Pole, a ~ 10 % 
variation is seen even at 3.2 GeV in the and Pg fluxes for vertically down going directions. 
The fluxes of the and also show variations, but smaller than those of the Ue and Pg. 
This is considered to be due to the change of muon energy loss rate in the atmosphere 
according to the change of the air density by the seasons. 

We also studied the production height of atmospheric neutrinos with the NRLMSISE- 
00 atmospheric model. For the mid-latitude region (SK site), they are very close to the 
production height calculated with the US-standard ’76 atmospheric model. In the polar 
region, we And a seasonal variation of production height, but in the mid-latitude and tropical 
regions, the seasonal variation of production height is small, and almost invisible. However, 
we And a large azimuthal variation of production height due to the muon bending in the 
mid-latitude and tropical regions. 

Here, we would like to make a short comment on the calculation error of the atmospheric 
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FIG. 15. Production height of atmospheric neutrinos for and I'e + going to vertically 

down and horizontal directions, summing all azimuth directions at the Pyhasahni mine. The height 
that cumulative distribution reaches 10 %, 50 %, and 90 % are shown as a function of neutrino 
energy. The time average values in June - August are shown as solid lines and in December - 
February as Dashed line. 


neutrino flux, and on the recent observations of the cosmic rays. First note that this work is 
still within the calculation scheme established in Ref which is based on the comparison 
of atmospheric muon observation and calculation. Therefore, the estimation made in Ref 
is still valid in this work. The total error is a little lower than 10 % in the energy region 
1-10 GeV. The error increases outside of this energy region due to the small number of 
available muon observation data at the lower energies, and due to the uncertainty of Kaon 
production at higher energies. 


The recent observations of the primary cosmic rays 17 h 20|. suggest that we need to 


modify the interaction model again to reconstruct the observed atmospheric muon fluxes 
with them. In a preliminary work 2]| , we hnd it is possible to modify the interaction model 
in such a way. and the atmospheric neutrino flux calculated with that is ve^ close to the 
present calculation in 1-10 GeV and well within the error estimated in Ref j5|. Now would 
be the time to study the interaction model with the updated primary cosmic ray spectra 
model for the calculation of atmospheric neutrino flux. 
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FIG. 16. Azimuthal variation of production height of atmospheric neutrinos at the SK site, INO 
site, South Pole and Pyhasalmi mine for near horizontal directions. The height that cumulative 
distribution reaches 50 % is shown as the function of azimuth for neutrino in the energy range of 
3.16 <E^< 3.55 GeV. 
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